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(<— ^ , The performance of the LAND neutron detector is studied. Using an event-mixing technique 

PsJ ' based on one-neutron data obtained in the SI 07 experiment at the GSI laboratory, we test the 

^-H \ efficiency of various analytic tools used to determine the multiplicity and kinematic properties 

^ ' of detected neutrons. A new algorithm developed recently for recognizing neutron showers from 

• '~j \ spectator decays in the ALADIN experiment S254 is described in detail. Its performance is 

rN ' assessed in comparison with other methods. The properties of the observed neutron events are 

j^ \ used to estimate the detection efficiency of LAND in this experiment. 
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1. Introduction 

The Large Area Neutron Detector (LAND) was constructed and used at the GSI Laboratory 
in Darmstadt for the measurement of neutrons from reactions at relativistic energies 111]. It con- 
sists of 200 modules (paddles) of 2 m length and 10x10 cm^ cross section, constructed from 
interleaved 5 mm thick iron and plastic-scintillator plates. They are arranged in 10 consecutive 
planes of 20 paddles each, alternating between horizontal and vertical paddle orientation, to form 
a compact 2x2x1 m"* detector 

The interaction of a neutron with the detector material results in the production of light pulses 
(hits) in one or more paddles. Photomultipliers register the light at both paddle ends, providing 
information on the amplitude and time of occurrence of each hit. The hit position along the 
paddle is determined from the time difference of the two signals. Each hit generated by a neutron 
is thus characterized by its position in space and time. Typical internal resolutions of the position 
and time measurements are approximately +3 cm, and ±0.25 ns, respectively ||2|]. 

The measured amplitudes are summed into the so-called visible energy, i.e. the energy de- 
posited in the scintillator material of the detector It is only a small part of the total energy 
deposited by an interacting neutron, as the majority of it is lost in the iron layers, invisible for 
the phototubes. Charged particles impinging on the detector are recognized by their signals in a 
VETO wall of 5-mm thick plastic-scintillator slabs mounted in front of the detector [[IJ]. 

LAND was mainly intended to serve as a time-of-flight spectrometer with the precise deter- 
mination of the position and arrival time of identified neutrons being the primary goal. Because 
of its calorimetric properties, it has also been used to study multi-neutron emissions in specific 
experiments. For determining the multiplicity of detected neutrons, two calorimetric observables 
can be used: the hit multiplicity A^hit and the total visible energy fivis, provided the average values 
of these quantities per incident neutron are known yj]. To obtain also the kinematic properties of 
individual neutrons, a more complicated data treatment is needed. It requires the identification 
of the first hit generated by the neutron, called here primary hit. Its registered location and time 
can be used to determine the velocity vector of the incident neutron. For neutron (or, more gen- 
erally, particle) multiplicities exceeding one, the individual hit patterns may overlap, turning the 
identification of primary hits into a complex problem. In fact, if the overlapping hit patterns are 
also generated close in time, unique solutions may not always be possible. A procedure for an 
optimized primary-hit recognition on an event-by-event basis is thus a necessary tool for experi- 
ments requiring multi-neutron detection, but also an interesting challenge. The assessment of its 
performance under the given experimental conditions will be essential. 

The standard neutron recognition algorithm used so far [J4|| seems to be successful in exper- 
iments where the neutron multiplicity to be expected is low. They include structure studies of 
halo nuclei as, e.g., "Li by fragmentation at 280 MeV/nucleon ||5|] or the excitation of collective 
modes as, e.g., the two-phonon giant dipole resonance in heavy nuclei la] or the pygmy dipole 
resonance in neutron-rich nuclei [|7|l. The same algorithm was also used for analyzing the data 
from the collective-flow experiments providing evidence for the neutron squeeze-out in '^^Au + 
'^^Au collisions at 400 MeV/nucleon ISllst] even though the most probable neutron multiplicity 
was of the order of 3 to 4, with nearly the same number of light charged particles being present 
in the recorded events. Because of considerably larger neutron multiplicities, in the ^^'^U + Fe at 
1 GeV/nucleon experiment, the neutron yields were measured with the use of only two forward 
LAND planes, with a very limited neutron detection efficiency |J2[]. 

Since similarly large neutron multiplicities were expected in the more recent study of projec- 
tile fragmentation at 600 MeV/nucleon, performed with beams of '"^''^"^Sn, '^'^La on ""'Sn targets 
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lIOllllllj a new algorithm for neutron-shower recognition was developed. It exploits the hit-to- 



hit time and space correlations expected for a shower and is intended to be especially adapted to 
the high neutron multiplicities from the de-excitation of highly excited projectile residues. The 
algorithm has been already used for the LAND data analysis in the S248@GSI and S304@GSI 



experiments performed by the SPALADIN Collaboration Ill2[ 11311 . It will also be available for 
the data collected in the very recent study of the density dependence of the symmetry energy by 
comparing the collective flows of neutrons and light charged particles from heavy-ion collisions 
at 400 MeV/nucleon. In this experiment of the ASYEOS Collaboration, the LAND detector 
served as the main tool for the measurement of the emissions of neutrons and hydrogen isotopes 
near mid-rapidity lll4ll5ll . 

Considering the wide range of applications and the complexity of the problem, it is indispens- 
able to acquire a quantitative knowledge of the performance of the used algorithm under the con- 
ditions of the particular experiment. Since ambiguous hit patterns without unique solutions have 
to be expected, this is even more important and of general interest for users of LAND. The present 
work is intended to serve this purpose by performing a comparative test of the standard and new 
algorithms with synthetic events generated from single neutron data from deuteron breakup. This 
experimental data set had been collected 20 years ago by measuring the response of LAND with 
tagged neutrons generated by dissociating deuterons of energies 70 to 1050 MeV/nucleon and 
by recording the protons deflected by the ALADIN magnet (experiment S107@GSI, see ref. 
llallfil])- Earlier tests with the same technique of generating synthetic multi-neutron events had 
shown that the quality of neutron detection and kinematic identification changes with the exper- 
imental conditions JlTIl . It is clear that it will also depend on the actual status of readiness of the 
detector; missing channels, e.g., will deteriorate the performance. Therefore, to make the simu- 
lation as realistic as possible, the measured hit distribution in time-space obtained in the '^"*Sn -i- 
Sn reaction at Slab = 600 MeV/nucleon, and the state of readiness of LAND in the frame of 
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the S254@GSI experiment, 1118142111 are taken into account. 

Even though this work is especially devoted to the LAND detector, the methods developed 
here can also be considered in the wider context of large neutron detectors. They should be 
useful for the NeuLAND detector proposed within the frame of the R3B project at the future 



FAIR facihty t2M 



2. Neutron recognition procedures 

2.1. Shower Volume Algorithm (SVA) 

The standard algorithm for neutron recognition in LAND, Shower Volume Algorithm, has 
been already described in 141 |5l [iTI] . Here we present a short description of the method, to give 
the reader an idea about the differences between the two algorithms studied in the paper 

In SVA the main focus is put on the time sequence of registered hits. The first registered hit 
is a primary hit, defining properties of the first detected neutron. Subsequently registered hits are 
attributed to the same neutron if they are located in a cylindrical volume of 29 cm radius and 
50 cm depth around the primary hit, or at a larger distance, if they fulfill the kinematic conditions 
of quasi-free scattering. The hits attributed to the neutron in this way are removed from further 
analysis. Then the same treatment is applied to the remaining hits. This will select a number of 
primary hits which, in the ideal case, will correspond to the neutrons that have hit the detector 

The cylindrical volume used in this procedure is quite large. Although it amounts to only 
3 % of the total LAND volume, one should remember that in a real experiment the neutron 



angular distribution may be strongly anisotropic so that only a small part of LAND is irradiated 
by the majority of neutrons. For this reason, the SVA algorithm starts to fail when the neutron 
multiplicity increases: the cylindrical volume frequently contains hits from more than only one 
neutron shower As a consequence, the resulting neutron multiplicity becomes underestimated. 
Moreover, as the hit-time order plays an important role, the faster neutrons are favored by the 
algorithm. In the events, where the neutron multiplicity is considerably large, this fact may 
disturb the observed neutron velocity distribution, underestimating the yields in the range of 
lower velocities. These disadvantages of the algorithm were the reasons for searching for a more 
refined algorithm devoted especially to the reactions characterized by a large neutron multiplicity. 

2.2. Shower Tracking Algorithm (STA) 

As mentioned above, a high-energy neutron or a charged particle, interacting with the LAND 
detector, produces usually several hits. To correctly determine the kinematic properties of the 
incoming particle one has to recognize properly the first ("primary") hit. Its time and position 
determine the velocity of the particle. The consecutive secondary hits are correlated in time and 
space with the primary hit, or a previously identified secondary hit, forming a kind of "jet" or 
"chain" of hits. Such a secondary hit identification can be done using recursive programming 
methods. This is the main idea of the Shower Tracking Algorithm. 

Let us consider two hits observed in two neighboring LAND planes, with time-space coordi- 
nates (xi,yi,ti) and ix2,y2,t2)- If the two hits are produced by the same interacting particle, they 
should be relatively close to each other Indeed, this is clearly observed when one studies plots 
of3'2 -yi vs. X2 -xi (Fig.[TJ. 

Analyzing hit-hit correlations one can define the correlation condition: 

-19 cm < Ax < 19 cm, 

-19 cm < A}? < I9cm, (1) 

-1 ns < Af < 3 ns. 

Using this condition one can identify a "cluster" of hits generated by a particle. A recursive 
procedure starts from a hit registered in the first LAND plane . For this hit, a search is made for 
a correlated hit in the same and in the next plane. If a correlated hit is found, a search for a hit 
correlated with this secondary hit is made. The procedure is continued until no more correlated 
hits are found. All correlated hits are marked as secondary hits and removed from the total pool 
of hits. Then the procedure starts from the beginning with one of the remaining hits, descending 
from the first plane to the second plane, third plane, and so on. 

After the clustering procedure has been completed, a number of "single" hits remain, where 
no correlation with other hits was found. Some of them are really "singlets", due to a small 
momentum transfer between the reacting neutron and nuclei, or when the reaction occurs close 
to the edge of the LAND detector But a number of them is a result of detection problems. In 
this number one should mention: 

(i) Missing paddles - a number of LAND paddles are not working (or improperly working) in 

the experiment. In this case some hit chains are broken. 
(ii) Multi-hits - when two hits are registered in one paddle, they are seen as one hit with a false 

position. Also its time is earlier than the real times, 
(iii) The neutrons can interact with LAND in two points, distant from each other These two 
points may be seen by the procedure as two separated clusters but should be attributed to 
one neutron. 
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We are not able to take properly into account the points (jl) and (jn)). We can only estimate their 
influence on the neutron multiplicity. To do this, we construct two limiting procedures. The first 
one, called "MAX", maximizes the neutron multiplicity: it treats all the remaining single hits as 
separate clusters. The second one, called "MIN", minimizes the multiplicity of neutrons, trying 
to clusterize the remaining single hits using a looser condition - only a wide time correlation is 
here taken into account: 

-7ns<Af<4ns. (2) 

A third procedure, called "AVE", is proposed to get a probably more realistic value in between 
these extremes. Here, for the remaining single hits, a condition is used that is more restrictive 
than ^ but still less restrictive than ([TJ: 

-19 cm < (Ax or Ay) < 19 cm, 
-1 ns < Af < 3 ns, 

i.e. only the x or the y separation is required to meet the correlation criterion. In the "MIN" and 
"AVE" procedures, the conditions (12} and ([3}, respectively, are applied in place of ([TJ in the main 
clustering procedure, and a second run of the procedure is performed for the ensemble of single 
hits. 

To take into account the point (Iml i we use a "long-range correlation procedure" searching for 
distant clusters correlated in time and space. For each previously determined primary hit we use 
its measured time and position to estimate the velocity of the incoming neutron. Then we search 
for another primary hit in a cylinder defined by this velocity vector as the cylinder axis, and a 
fixed radius of 19 cm. Such a hit is treated as correlated with the first hit, if its time fulfills the 

relation: 

d 



Ar-- 

V 



< Ins, (4) 



where d is the distance between the two primary hits and v is the velocity of the first primary hit. 
This treatment is added as the last step in the procedures "AVE" and "MIN". 

In the following, the STA acronym refers to the "AVE" procedure, which will be shown to be 
the most accurate. The "MIN" and "MAX" procedures are considered as error estimators. 

3. Testing efficiency of neutron recognition 

3.1. Mixing single-neutron data 

The validation of the neutron recognition procedures has been performed by analyzing syn- 
thetic multi-neutron events generated from single-neutron data provided by the S107@GSI ex- 



periment III6I1 . In the experiment the deuteron beam was scattered on a ^'*^Pb target at energies 
from 70 to 1050 MeV/nucleon. The experimental set-up was focused on deuteron break-up 
registration. The neutrons hitting LAND were measured in coincidence with the accompanying 
proton, deflected by the magnetic field of the ALADIN magnet. In this way the information 
about the showers produced by a single neutron in LAND was collected. From this experiment, 
an efficiency curve for neutron detection with LAND was obtained 161 [iTI] . For neutrons of 600 
MeV, e.g., a detection efficiency of 94% was observed with the fully operational detector More- 
over, the single-neutron events registered in LAND can be randomly mixed at random positions, 
simulating in this way two-, three-, or more-neutron events with predetermined time and spatial 
distributions. 



We used this technique to determine the accuracy of the neutron recognition algorithms and 
the systematic and statistical errors generated by various analytic tools. In the mixing procedure 
we take into account the treatment used before in yl] for multiple hits in one paddle: 

1 . A new time is calculated by taking into account only the signals arriving first at the two 
photo-tubes, and then the new hit time and position are derived (this is a slight approxima- 
tion to the actual response of the electronic readout when the signals overlap); 

2. The signal energies of the two hits are summed at both paddle ends, and then the visible 
energy is recalculated. 

Furthermore, the actual spatial and temporal hit distributions, as measured in S254@GSI 
experiment, are reproduced in the mixed data. Also non-working paddles in this experiment 
have been removed from simulated events. 

3.2. Resolving efficiency and quality 

In the evaluation of a recognition algorithm one should consider two main problems: some 
neutron showers detected in LAND are lost due to low time-space resolution of the applied algo- 
rithm - i.e. two or more showers are erroneously recognized as a single one. On the other hand, 
in some cases the algorithm produces "false" showers, i.e. a one-neutron shower is artificially 
identified as two or more showers - because of either inadequacies of the algorithm, or some 
detection problems. From the point of view of neutron recognition, only the assignment of the 
primary hit is important. A mistake in attributing secondary hits to appropriate showers does not 
play any role. 

Besides the correct determination of the actual multiplicities, the recognition of the individual 
neutron showers is also important. A minimum requirement is that all hits, identified as primary, 
are actually belonging to different showers. In the simulated multiple neutron data it is possible 
to check the origin of each primary hit recognized by an algorithm. If a number n of them belong 
in fact to the same neutron shower (i.e. the same S107@GSI event) we speak about one "true" 
neutron and n - 1 "false" neutrons. If a neutron shower does not have a representative among the 
primary hits, we call this neutron "lost". 

Let us denote the number of true neutrons as Nt, and the number of lost neutrons as Nl- 
Then the total number of initial neutrons A^ can be expressed by the sum: 

Nt+Nl^ N. (5) 

On the other hand, the total number of recognized primary hits Nc can be divided into two 
groups: the number of true neutrons Nt, and the number of false neutrons Nf: 



Nt+Nf^ Nc- (6) 



Let us define the variables: 



Nt 
A = — , (7) 

N 

Nt 
Q = ,, ,, ■ (8) 

Nf +Nl 

The resolving efficiency A expresses the rate of correctly identified neutrons. It is evident from 
Fig. 12^ that, independent of A^, STA is more efficient in recognizing true showers than SVA. But 
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this is not a sufficient argument to favor STA. Let us imagine an algorithm treating each registered 
hit as a separated neutron shower Its resolving efficiency will be unity for all A^, because A does 
not carry the information about the number of false showers Nf- A more useful variable for 
characterizing an algorithm is the resolving quality Q. The larger its value is, the better the 
algorithm works. Fig. |2j) shows the dependence of this quantity on neutron multiplicity. At low 
neutron multiplicities, the resolving quality of SVA is higher but it decreases dramatically when 
the multiplicity becomes larger Above the crossing point at A^ = 6 the resolving quality of STA 
is superior For this reason, STA should be applied in experiments with high neutron multiplicity. 
When the expected neutron multiplicity does not exceed 6, SVA appears to be the better choice. 

3.3. Neutron multiplicity estimation 

We consider here four estimators for the determination of the neutron multiplicity on an 
event-by-event basis: 

(a) Normalized hit multiplicity A^hit//^, 

(b) Normalized total visible energy fivis/f, 

(c) Primary hit multiplicity given by SVA, 

(d) Primary hit multiplicity given by STA. 

In observables daj and ©, fi and e are, respectively, the mean hit multiplicity and the mean visible 
energy generated by one incoming neutron. These parameters can be determined from one- 
neutron data. Unfortunately, the results of the S107@GSI experiment cannot be directly used in 
other experiments, as specific setup conditions (in particular CFD thresholds and geometry) can 
considerably change these mean values. Alternatively, one can apply SVA and/or STA algorithms 
to select one-neutron events from the studied experimental data. We apply this latter concept to 
determine the values fi and e and we adopt the average of the two values given by the SVA and 
STA algorithms. Such a treatment used with the S107@GSI data has given the results which 
differed only by few percent from the actual ju and e derived directly from the data. 

The estimated neutron multiplicities are compared with the known multiplicities of synthetic 
multi-neutron events in Fig. |3] The solid lines represent the average real multiplicity as a func- 
tion of the estimated value. The distribution of the symbols reflects the uncertainties of the 
procedures. The evidently most accurate result is obtained with the estimator based on the ^vis 
observable. This is a simple consequence of the fact that in the mixing of one-neutron events we 
simulate the multi-hit visible energy by summing all partial energies. If this model is correct, the 
total visible energy is retained in the case of double hits, whereas the number of hits does not 
grow if a paddle sees two hits. The resulting upbend of the true vs. estimated multiplicity can be 
calibrated but the resolution becomes poor at high multiplicity (Fig. [3^). It is also evident that 
the proportionality factor e is quite well determined by the recognition algorithms. 

In the lower row of Fig. |3]one can compare the performance of the two recognition algo- 
rithms. It is evident that SVA fails in neutron multiplicity determination when the multiplicity is 
larger than 4. 

A more precise study of the errors generated by the estimators is presented in Fig. |4] The 
relative systematic error (RSE) and relative standard deviation (RSD) can be expressed as follows: 

RSE = ^^, (9) 

RSD = --!—^ . (10) 
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Fig. |4^ shows the comparison of RSE as a function of the estimated neutron multiplicity. It can 
be seen that the calorimetric estimators give relatively better accuracy for larger multiplicities. 
A similar result can be observed in Fig. |4j) where RSD is presented. In the range of lower mul- 
tiplicities, the recognition algorithms give better accuracy in comparison with the calorimetric 
estimators. For larger multiplicities fiyis gives the best estimate. 

Concluding, when only the total neutron multiplicity is needed, the fivis/e estimator gives the 
best measure. A typical systematic error in an event-by-event analysis is lower than ±10 %. It 
corresponds to an absolute neutron multiplicity deviation better than + 1 . The statistical errors 
are below 25 % for multiplicity 8 and larger To improve the predictions in the range of small 
multiplicities one can combine this estimator with the others. 

Let us study the predictive power of the estimators for an entire multiplicity distribution. 
Using mixed neutron data one can simulate a Poissonian multiplicity distribution and compare it 
to the distribution reconstructed with different multiplicity estimators. Fig.|5]shows the results of 
such a simulation done for the neutrons from the S254@GSI experiment at 600 MeV/nucleon. 
Evidently, SVA reproduces the multiplicity distribution very well if the mean multiplicity is 
around 2. When the multiplicity increases, this algorithm is less efficient, and the use of STA 
will be more appropriate. The best estimate is again provided by fivis/f- 

In an actual experiment, a part of the LAND paddles may not be working. This was the 
situation in the S254@GSI experiment, where about 15 % of LAND paddles were out of order, 
mainly in the rear part of the detector The right column of Fig. |5] shows that the obtained 
multiplicity distributions are very similar, even in the case of such a defect. 

4. Neutron detection efficiency 

The efficiency of LAND for neutron detection is high for high neutron energies. For neu- 
trons from the deuteron breakup at 600 MeV/nucleon, a value 770 = 0.94 was measured in the 
S107@GSI experiment with the then fully functional detector In a given experiment, the actual 
detection efficiency may deviate from this value because of a specific geometry, the aging of the 
scintillator materiaU and possibly modified thresholds, or non-working paddles. The efficiency 
is also expected to change near the edges of the detector In the following subsections, these 
aspects will be addressed. 

4.1. High efficiency area 

The S254@GSI experiment had the layout of instrumentation nearly identical to that pre- 



sented in Fig. 1 of Ref. II23I1 . The LAND detector was positioned approximately 10 m down- 
stream from the target, vertically symmetric with respect to the beam axis but displaced hor- 
izontally by 78 cm. It resulted in a symmetric vertical {y direction) but strongly asymmetric 
horizontal (x direction) distribution of hit positions. The magnetic field of ALADIN does not 
affect the neutron trajectories, and the neutron emission is expected to have a cylindrical sym- 
metry, with similar hit distributions in both vertical and horizontal directions. Because of the 



finite acceptance of the ALADIN magnet II23I1 . shadowing effects near the edges of LAND are 



expected. Examples of the horizontal and vertical hit distributions obtained from the horizontally 



'The reduction of the average attenuation length from originally 215 cm Q|| to 116 cm at the time of the S254 
experiment, is evidence for a degradation of the detector material with time. 
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and vertically oriented paddles, respectively, are shown in Fig. |6]in the beam-centered coordi- 
nate frame. The expected azimuthal symmetry is observed up to a distance of about 70 cm from 
the center. The more rapid vertical decrease beyond 80 cm coincides with the shadow of the 
ALADIN poles and magnet chamber as projected from the target position. The horizontal dis- 
tribution extends further out. Here the projected shadow of the magnet is close to the end of the 
detector Both distributions start to drop more rapidly at about 15 cm from the end of the detector 
but extend beyond it because of the finite position resolution obtained from the time differences 
at the paddle ends. 

Also the yield observed on the opposite end of the horizontal paddles (at -)-20 cm) is lower 
than that observed in the corresponding section of vertical ones. It suggests some additional edge 
effect, maybe connected to lower detection efficiency at the ends of a paddle. 

Because of these obvious edge effects and in order to work with a homogeneous detection 
efficiency we select a high efficiency area by cutting off the edges of the detector. This area is 
defined by: 

-150 cm < .Y < cm, 

-70cm <y< 70cm. ^ ' 

4.2. Estimation of missed neutron fraction 

The distribution of the number of hits Fo(^), generated by single neutrons of energy 600 Me V 
from the S107@GSI data and for the fully operational detector, is shown in Fig. |7]by the solid 
histogram. The measured efficiency for this energy is ^o = 0.94 16|], indicating that 6% of the 
incident neutrons pass the detector without a visible interaction. The corresponding probability 
yo(0) = 0.06 is given by the dashed line in Fig. |7] 

The hit multiplicity distribution for A^ = 1 in the S254@GSI experiment, Yf.^p{k), obtained by 
selection of single neutron events with the procedure "AVE" (with primary hit located in the high 
efficiency area) is represented by the points in Fig. |8] Comparing with yo(^), the distribution 
is much narrower, the average number of hits per neutron is almost halved. This significant 
reduction of the hit multiplicity is mainly due to non-working segments of the detector (~ 15%) 
and a considerably higher hit registration threshold. 

The estimation of the detector efficiency requires an extrapolation of the hit distributions to 
include their value for ^ = 0, assuming a realistic functional form of the distribution. 

First, let us apply the Poissonian distribution ypois(^) - Cj^e^'^ with two adjustable param- 
eters: the normalization constant C and the mean value A. The parameters are determined from 
the conditions Ypois{k) - Ye^j,{k) for ^ = 1 and k - 2, that provide the best approximation in 
the low multiplicity range. The resulting Fpois(^) distribution with A — 1.02 is shown in Fig. [8^ 
by the histogram. With the assumption yexp(O) = ypois(O) the procedure suggests the efficiency 
?7= 1.0 -exp(-1.02) = 0.64. 

Alternatively, one may try to determine Y{k) with a model closer to the actual detector sit- 
uation. We use the distribution Y^ik) measured with the fully operational detector as a primary 
distribution which is then modified by rejecting a certain fraction of hits. 

The hits are considered as a linear chain. Each hit in the chain is accepted or rejected with a 
given rejection probability P. If the chain is broken by rejected hit(s) into two or more separate 
parts, such an event is removed from the considered sample of events with N - I (the clustering 
procedure would recognize two or more neutrons). After applying this random rejection method 



the final distribution is found as: 

Yfik^O) = CfYjP'Yod), (12) 

i>0 

Yf(k>i) = Cf(i-p/2o- + i-^)P'-*yo(0. (13) 

i>k 

Requiring Y{(k) - Yf^^p{k) for ^ = 1 and k = 2 gives Cf and P = 0.507 with the result shown in 
Fig- IS'- Using the expression for k - (Eq. (fTSI) ). the efficiency is found to be 77 = 0.73 for this 
particular experiment. This is significantly lower than the value of 94 % previously deduced from 
the deuteron-breakup experiment S107@GSI performed with the then new and fully functional 
detector The estimated detection efficiency is specific for the S254@GSI experiment and may be 
different in other experiments. With the presented method, it can be obtained from the measured 
data. 

5. Conclusions 

An overview has been presented on the analytic tools, applied so far in the analysis of the 
LAND data. We compare two calorimetric observables and two different neutron recognition 
algorithms. The work estimates the systematic and statistical errors for an event-by-event neutron 
multiplicity estimation. The new algorithm, recently developed for the S254@GSI experiment, is 
described in detail and compared to the algorithm used before. A discussion of the applicability 
of each algorithm is presented. The new algorithm seems to be a good choice for the analysis 
of the data, when expecting neutron multiplicities larger than 4, especially when the emission of 
both neutrons and protons has to be studied. 

The study allowed to estimate the neutron-detection efficiency of the LAND to be about 
70 %. It is worth noting that this estimate is significantly lower than the value 94 % deduced 
many years ago from the results of the deuteron-breakup experiment S107@GSI. The difference 
is believed to be due to the natural aging of the scintillator material, possibly higher effective 
thresholds of the discriminators, and to the about 15 % non-working paddles. This emphasizes 
the need for methods for determining the detection efficiency from the actual data measured in a 
specific experiment. Significant variations can, apparently, not be excluded and will be difficult 
to model since some of the required parameters may not be sufficiently well known. The methods 
presented here can be expected to serve this purpose. 
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Figure 1: A typical iiit-iiit correlation for two hits registered in two neighboring planes and in the time interval -1 ns < 
r2 — ?i < 3 ns (experimental data from S254@GSI) 
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Figure 2: Resolving efficiency (a) and quality (b) obtained for STA (closed points) and SVA (open points) algoiithms, as 
functions of the real neutron multiplicity. The missing neutron fraction (see Section |4^ is not taken into account. The 
simulation uses actual spatial and temporal hit distributions measured at £lab = 600 MeV/nucleon incident energy in 
the S254@GSI experiment 
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Figure 3: Comparison of neutron multiplicity estimators: (a) A'hit/jU, (b) i?vis/f. (c) SVA, (d) STA. The solid lines show the 
average neutron multiplicity as a function of the estimated value. The boxes represent the event distributions. The gray 
lines in (d) correspond to the "MIN" and "MAX" limiting procedures. The thin solid line shows the 1:1 coiTespondence 
case. The simulation uses actual spatial and temporal hit distributions measured at £lab = 600 MeV/nucleon incident 
energy in the S254@GSI experiment 
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Figure 4: Comparison of relative systematic errors RSE (a) and relative statistical deviations RSD (b) generated by 
neutron multiplicity estimators, as functions of the estimated neutron multiplicity. The simulation uses actual spatial and 
temporal hit distributions measured at £lab = 600 MeV/nucleon incident energy in the S254@GSI experiment 
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Figure 5: Simulated Poissonian neutron distributions (circles) and reconstructed distributions using £vis/f estimator 
(solid line), STA (dashed line) and SVA (dotted line). The left column corresponds to full LAND performance, the right 
column to only 85% working paddles. The mean neutron multiplicities ai'e: N = 2 (upper row). N = 6 (middle), and N = 
10 (bottom). The simulation uses actual spatial and temporal hit distributions measured at Slab = 600 MeV/nucleon 
incident energy in the S254@GS1 experiment 
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Figure 6: Distribution of absolute hit positions along a paddle for horizontal (black line) and vertical (gray line) paddles, 
as observed in the S254@GSI experiment. Dashed vertical lines show the corresponding geometrical limits of the 
detector. Triangles indicate the software limits l lll) imposed on the data to get a homogeneous efiiciency ai'ea. The 
recorded hit-position distributions exceed the actual dimensions of LAND because of the finite resolution of the time 
measurement 
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Figure 7: The hit multiplicity distribution for the fully operational detector (solid line) and estimated fraction of missed 
neutrons (dashed line). Result obtained from S107@GSI experiment 
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Figure 8: The experimental hit multiplicity distribution for neutron multiplicity N = \ with its approximations by a 
Poissonian distribution (a) and a distribution from the hit rejection procedure (b). The lower panels show the same in log 
scale 
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